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The proximity of antiferromagnetic order in high-temperature superconducting materials is con- 
sidered a possible clue to the electronic excitations which form superconducting pairs. Here we 
study the transverse and longitudinal spin excitation spectrum in a one-band model in the pure 
spin density wave (SDW) state and in the coexistence state of SDW and the superconductivity. We 
start from a Stoner insulator and study the evolution of the spectrum with doping, including dis- 
tinct situations with only hole pockets, with only electron pockets and with pockets of both types. 
In addition to the usual spin-wave modes, in the partially gapped cases we find significant weight 
of low-energy particle-hole excitations. We discuss the implications of our findings for neutron 
scattering experiments and for theories of Cooper-pairing in the metallic SDW state. 
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I. INTRODUCTION 

Understanding the microscopic origin of interacting and competing phases with multiple order parameters (OP) is 
one of the central issues in condensed matter. This topic is particularly important for describing the complex phases 
of layered transition metal oxides such as the high-T c superconductors in which superconductivity occurs upon hole 
or electron doping of an antiferromagnetic Mott insulator. As a consequence of the competition between the Mott 
insulator and superconductivity, a rich variety of ground states may emerge, such as charge, spin, or current density 
wave states, all of which have been proposed to play a role in the pseudogap region of the high-T c phase diagram. 
For n-type cuprates, there is growing experimental evidence for the evolution of superconductivity in the background 
of the commensurate spin density wave (AF) state with ordering momentum Q= (jr, 

One of the most interesting questions concerns the evolution of the spin excitations in a pure AF state upon 
introducing carriers by doping, and the possible coexistence of such a state with d x i _ y i -wave superconductivity. The 
former question was studied in detail within a strong-coupling models such as the t — J model Hamiltonian^H It starts 
from a model of localized spins at zero doping which interact via the antiferromagnetic superexchange interaction 
and then introduces doped holes on the AF background. This approach is well justified in the hole-doped cuprates 
which are Mott insulators in the AF state. At the same time, experimental studies of the moderately electron-doped 
cuprates 1 and parent iron-based superconductors 6 show that both remain metallic in the AF state. In addition, a 
recent DMFT study of the undoped electron-doped cuprates with T' crystal structure 7 suggests that their insulating 
behavior is more due to the presence of magnetic long-range order than to charge transfer gap physics. This implies 
that a conventional weak-coupling description of the spin excitations may be applied for the n-doped case. This 
approach was initially developed in the nine ties^ , and more recently shown to describe well the spin wave spectra of 
the parent cuprate superconductor La2CuOi^^ and iron-based superconductors^. At the same time, less is known 
about the spin excitation spectra in the doped situation when the AF state in the electron-doped cuprates is still 
metallic. 

In this paper, we employ the RPA formalism within a single-band Hubbard framework to describe the evolution of 
the spin excitations from a Stoner insulator to the antiferromagnetic metal. In particular, we analyze the evolution 
of the transverse and longitudinal spin excitations when the system experiences a transition from (i) the Stoner 
insulator to the antiferromagnetic metal with electron Fermi surface pockets centered at (±7r, 0) and (0, ±7r) points 
of the Brillouin Zone, and (ii) the Lifshitz transition between two types of the AF metal when in addition to the 
electron pockets the second small hole Fermi surface pockets appear at (±7r/2, ±7r/2) points of the Brillouin zone (BZ). 
Furthermore, motivated by the observed coexistence of the AF and d x 2_ y 2-w&ve superconductivity in the electron- 
doped cuprates we compute the spin excitations in the coexistence region. We find that the Goldstone mode in the 
transverse channel remains gapless in the coexistence regime and that the excitations in the transverse channel are 
dominated by the renormalized spectrum of the spin waves. At the same time, we find that the excitations in the 
longitudinal channel include a resonance mode at the commensurate momentum close to (ir, 7r). The simultaneous 
coexistence of the longitudinal resonance and the transverse spin waves opens up an interesting possibility to use 
inelastic neutron scattering (INS) to identify the microscopic coexistence of superconductivity and antiferromagnetism. 
We also stress the importance of well-defined particle-hole branches of the spin excitation spectrum which should also 
be observable in INS and should help to confirm our general picture. 
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FIG. 1: color online) Three possible types of Fermi surface topology in the commensurate SDW state in layered cuprates either 
for hole (a) or electron doping (b) and (c) analyzed in this paper. Due to the mean-field SDW order the original large Fermi 
surface (green curve) folds, yielding small hole pockets (blue curves) centered around (±7r/2, ±7r/2) and electron pockets (red 
curves) centered around (±7r, 0) [((0, ±71"))] points of the BZ. For larger doping and smaller sizes of the AF gap both types of the 
pockets can be present. As argued in the text, the commensurate SDW order becomes unstable once the hole pockets appear 
around (±7r/2, ±7r/2) points of the BZ. This occurs due to overall negative spin stiffness of the commensurate spin waves in 
this case, see Figpna). 



Our starting point for investigating the coexistence of superconductivity and SDW order is the Hamiltonian 
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where c£ (cka) creates (annihilates) an electron with spin a and momentum k. We consider a two-dimensional 
system with normal state tight-binding energy dispersion = — 2t (cos k x + cos k y ) +4t' cos k x cos ky—fi and hopping 
matrix elements between nearest (t) and next-nearest (t') neighbors. The chemical potential \i controls the doping, 
x which is determined by n = 1 + x. The second and third terms in Eq. lead in the appropriate Hartree-Fock 
factorization to the emergence of commensurate SDW order and superconductivity, respectively. While it is generally 
assumed that both phases emerge from the same underlying interaction, its renormalization due to vertex corrections 
gives rise to a different effective interaction in each channel, thus justifying the Hamiltonian of Eq.QP 

We first focus on the pure SDW ordered state, and consider the case Vq = in Eq. (lj. After decoupling the second 
term via a mean-field (MF) approximation and diagonalizing the resulting MF Hamiltonian via a conventional SDW 
Bogoliubov transformation, we obtain two electronic bands (labeled a and /3) in the reduced Brillouin zone (RBZ) 

with dispersions — ± \j (e^ ) 2 + W 2 where W = cr (cj c / + Q a c\ c ^ a )sgna is the SDW order parameter 

which is determined self-consistently for a given U, and = (£k±£k+Q) /2. For completeness, one also has to 
include the self-consistent determination of the chemical potential. Depending on the strength of the SDW gap and 
the doping level, the system can be either a (i) Stoner insulator for the half-filled case and large value of W; or (ii) 
SDW metal with small size electron and hole Fermi surface pockets for intermediate values of W and non-zero electron 
or hole doping, respectively. Typical Fermi surfaces in the SDW metal are shown in Fig. [l] 

The dynamical spin susceptibility for the longitudinal, zz, and the transverse, H — , components is defined as 
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where Im = zz,-\ — . As mentioned above, the antiferromagnetic ordering at Q= (7r,7r) doubles the unit cell and 
requires accounting for the breaking of translational symmetrySEE. As a result, the total susceptibility in the transverse 
channel is a 2 x 2 matrix 



x+-(q, q ,fi) x + "(q,q + Q,^) 
x + -(q + Q,q,^) x + "(q + Q,q + Q,tt) 



(3) 
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Then, within RPA the susceptibility is obtained via solving a Dyson equation and the result can be expressed^ as 
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and the bare components are given by 
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with 7 = a, /? and for the Umklapp term 
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where /(e) is the Fermi function and the prime refers to the sum over the magnetic (reduced) Brillouine Zone (MBZ). 
Observe that we have Xo" ~~ (q> Q + Q) = Xo" ~~ (Q + Q> Q) anc ^ the expression for Xo" ~~ (q + Q? Q + Q? ^) can be obtained 
from Eq. 

The Umklapp component is absent for the longitudinal susceptibility and the RPA expression is given by 



with the bare longitudinal susceptibility in the form 

/ 
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Here, we note the opposite sign of the W 2 term in the coherence factors of the longitudinal and transverse suscep- 
tibilities. Although the general structure of the spin susceptibility in the SDW state at Q is knowrP^, it is useful 
to mention its main features. In particular, below TV spin-rotational invariance is broken and x H ^X zz - The 
imaginary part of the transverse component is gapless and displays the Goldstone mode at Q = (7r,7r) and Q — >• 
for any temperature in the range < T < Tjy. The latter is guaranteed by the fact that the condition of the pole 
formation in the RPA part of the transverse spin susceptibility coincides with the mean-field equation for W and, 
as clearly seen, is valid for any doping level as soon as the SDW order exists. At the same time, the longitudinal 
component of the spin susceptibility at Q is gapped by twice the SDW gap magnitude, W. 
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FIG. 2: (color online) Calculated spin wave dispersion Q vs. q in units of tv/cl along the symmetry route (0,0) — >> (7r,7r) — >• 
(0,?r) (0,0) of the BZ for a Stoner insulator at zero doping (x=0), W = 0.75t and t' = (a), t' = 0.2t (b), and t' = 0.35£ 
(c). Here, we fixed the magnitude of W by using [7 = 1.40t and the chemical potential \i — O.OOt (a) \i — — 0.30t (b), and 
[i — — 0.68£ (c), accordingly. The intensity for Imx^p A (q, q, H) is shown on the log scale. 

What is less well known, however, is the behavior of spin excitations away from Q. At half- filling (x=0) and zero 
temperature, the Fermi surface is gapped for any value of W and t' enforced by the self-consistent calculation of the 
chemical potential.^ The excitations in the transverse channel of the Stoner insulator are spin waves - collective spin 
modes of the antiferromagnetic ground state - with a dispersion all over the BZ 10 . Due to full gapping of the Fermi 
surface, the particle- hole Stoner excitations and the spin waves are separated in energy and may interact only around 
Qp-z^q), whose magnitude is controlled by the values of the SDW gap magnitude W and t' jt ratio. For t' = the 
onset of the particle- hole continuum is gapped at least up to Cl p _h(Q) = 2W. This is because the top of the lower 
/3-band and the bottom of the upper a-band are located at the magnetic BZ boundary, i.e. cosk x + cosk y = 0, 
at energies — W and +W, respectively. Therefore, there exists a degenerate manifold of q wave vectors for which 
Q p -h(q) = 2W. As a result, the spin waves do not interact with the particle- hole continuum for sufficiently large 
values of W and look identical to those obtained within a Heisenberg model of localized spins which interact via an 
AF exchange between nearest neighbors, Ji ~ jj, see Figj^a). 

The situation changes for non-zero values of t' jt. Its existence (i) introduces non-degenerate positions of the top of 
the j3— band and bottom of the a-band, (ii) reduces the overall magnitude of the gap in the particle- hole continuum 
and shifts it to lower energies at the (7r/2,7r/2) point of the BZ as clearly seen from Figj2^b)-(c). To understand its 
origin, observe that for any non-zero t' jt the bottom of the upper a-band is located at (±7r, 0) and (0, ±7r) points of 
the BZ at energy —At' + W — n > 0, while the top of the lower j3— band is located at (±7r/2, ±7r/2) points of the BZ 
at energy — W — /i < 0. As a result, the smallest gap between both bands which determines also the lowest position 
of the particle-hole continuum occurs at Q p -h(q_i) = 2W — At' for qi = (±7r/2, ±7r/2). For increasing t' jt ratio and 
constant value of W, the spin waves are bounded from above at momentum qi = (±7r/2, ±7r/2) and form a local 
minimum there at energies below 2W — At' > 0. In particular, in Figj2|b) it occurs below 0.7t with W = 0.75t and 
t' = 0.2t and is shifted to much lower energies for t' = 0.35t for a fixed W = 0.75t as shown in Fig{2|c). Notice also 
that for zero doping we always find either an insulating SDW state (which guarantees the finite gap at momentum 
qi) or a paramagnetic metallic state for U < U c . 

The observation of the local minimum for a finite t' at qi is due to an interaction of spin waves with the particle-hole 
continuum. This is certainly a signature of weak-coupling as it requires 2W to be of the same order as At' . This 
effect would not occur for the localized J\ — J2 model, where J2 refers to the antiferromagnetic exchange between 
the next-nearest neighbors. The inclusion of J 2 usually only lowers the position of the maximum of the spin wave 
dispersion at the Y point of the BZ, an effect clearly reproduced in the weak-coupling calculations as well [compare 
Figj^a) and (c)]. At the same time, within the localized model the particle-hole excitations always remain gapped 
by the large value of U [W]. Correspondingly the local minimum in the spin susceptibility at qi never forms in this 
case. 

To finish the discussion of the spin waves for the undoped case, we further notice that the spin wave dispersion 
obtained is symmetric with respect to the (0, 0) and (jr, it) points, which reflects the fact that both are equivalent 
symmetry points of the magnetic BZ. At the same time, the absolute intensity of the spin waves is different and is 
determined by the SDW matrix elements which are suppressed around the T-point. Indeed, from Eq.([5| one finds that 
at low Q the non- vanishing contribution to the intensity comes from the interband (a — >> f3 and visa versa) transitions 



which are proportional to the SDW coherence factor c l ^ t ^ r 
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For q ~ Q one finds 
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FIG. 3: color online) Calculated transverse, Imx^ PA (q, q, H), (upper panel) and longitudinal, Imxfjp A (q, q, Q) (lower panel) 
spin excitation spectra Q vs. q in units of ir/a for the metallic commensurate SDW state with t' /t = 0.35 and U = 1.3875£. 
(a),(d) refer to the hole doping, x=-0.05, W = 0.61£, p = -0.8819*, (b),(e) refer to the electron doping , x = 0.10, W = 0.4404t, 
p — — 0.4284£, and (c),(f) refer to the electron doping x = 0.14, W = 0.12£, p — — 0.302£. The corresponding Fermi surface 
topology is shown in Fig. [I] The intensity in states /t is shown on a log scale. The white arrows in (a) denote the incommensurate 
momentum associated with 2k f scattering due to hole FS pockets shown in FigJTJa). Observe also the difference in the intensity 
maps between upper and lower panels. 

£ k+Q ~ - £ k and 4"q~Q ~ 2 > whereas it is cggZo oc (l - ( ( ^} 2 ~^ 2 ) for q - 0. 

Let us continue by looking at the spin excitations of the commensurate AF order upon electron or hole doping. As 
we are interested in the metallic ground state, we will only study the situation of the AF metal for finite t' jt ratio and 
remain at very low temperatures. We remind the reader that in this case the Fermi surface topology for the electron- 
and hole-doped AF metal is different. In the hole-doped case the chemical potential crosses first the top of the lower 
/3-band, which is located at (±7r/2, ±7r/2) points of the BZ, while for the electron doping the chemical potential shifts 
to the bottom of the upper a— band which results in the FS electron pockets at (±7r,0) and (0, ±7r). 

For the hole-doped case and t' jt > 0, we find that the commensurate antiferromagnetic order is unstable for any 
finite doping, consistent with previous analysed- 12 * 14 '. An instability of the commensurate SDW order can be clearly 
seen from Figj3|a) where we show show an example of the spin excitations obtained for the commensurate SDW order 
and x — —0.05. Note that despite a commensurate SDW order parameter the spin excitations form gapless dispersive 
incommensurate modes at in addition to the strongly damped Goldstone mode which exists at Q, despite the fact 
that commensurate ordering at (7r, 7r) was assumed in the mean field decoupling. This signals that the commensurate 
SDW order is not a stable solution. The origin of this instability is related to the appearance of the small FS hole 
pockets, whose presence introduces two effects. The first one is related to the spin stiffness of the commensurate 
spin excitations at Q. In particular, in contrast to the undoped case, there is an additional contribution to the spin 
stiffness which arises due to intraband ft — (3 transitions which are now gapless. To see this, we expand the dispersion 

2 2 

of the lower /?— band for U » t close to the (±7r/2, ±7r/2) points which yields E^ — —p — W — — where 

P\\ = (k x — k y )/2, p± = (k x + k y )/2 and m\\ = (St') -1 , and m±_ = (16t 2 /W — St') -1 . As follows from the analysis 
of the denominator of the transverse spin susceptibility at Q, the spin stiffness, p s acquires a finite correction in the 
doped SDW metaP^n3HU as p g = p®(\ — z) where z = 2U ^ /rn ^ m — is proportional to the Pauli susceptibility of the 
/3-band, and is the bare spin stiffness in the undoped case. As one clearly sees, z > 1 for large U which indicates 
that the commensurate SDW order is unstable upon hole doping. On the contrary, for the opposite case U « t the 

expansion yields E£ = —p — ^ v±P± + ^7 where here m± = m\\ = (S^) -1 and v± ~ t. This indicates that 

for t' < t the dispersion along p± is essentially linear. As a result the static susceptibility of the /3-band will have 
singular behavior at 2kp. Therefore, even if the Pauli susceptibility of f3 electrons yields z < 1, there is an additional 
source of instability associated with 2kp scattering on the hole pockets. To verify which of these two instabilities 
occurs for U ~ t, used in our numerical calculations, we analyzed the behavior of the denominator of the RPA spin 
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susceptibility and found that for the case of Figj3[a) we have z > 1. Therefore, the instability of the commensurate 
SDW order for U ~t and hole doping occurs due to negative corrections to the spin stiffness. 

For electron doping the situation is different. In this case the chemical potential moves towards the a— band which 
has a minimum at (±7r,0) and (0, ±7r) points of the BZ. Expanding now the dispersion of the a— band around them 

yields Eg = 4£' — M + ^+ ^ + ^m wnere m ~ (4t / ) _1 . Observe that this expansion holds for any ratio of U/t. As 
in the case of the hole pockets, there is a contribution of the a — a scattering to the spin stiffness. However, it is easy 
to see that for the intermediate values of U ~ t and t' = 0.35t one has z < 1. In addition, there are additional & 4 
corrections in the dispersions (proportional to t 2 /W) which further reduce the value of z. We verified it numerically 
by varying values of U and still finding commensurate SDW order to be a stable solution for electron doping. 

Thus, one of the main results of our paper is that for non-zero t' jt ratio there is a strong anisotropy in the evolution 
of the commensurate SDW order in the hole- and electron-doped cases [Note that the situation is just opposite 
for the negative sign of t'\. While for hole doping the commensurate SDW order is unstable due to an additional 
response associated with the hole pockets at (±7r/2, ±7r/2) which introduce negative spin stiffness to the spin waves, 
the influence of the electron pockets appearing at (±7r, 0) or (0, ±7r) is less dramatic and the commensurate SDW 
persists to a much larger doping. Observe, for example, that for x = 0.1 (Fig/j3jb)) the spin waves disperse all over 
the BZ despite the damping introduced by the a — a scattering which is clearly seen in Figj3^e). It is tempting to 
connect our findings to the experimental observations which show that the commensurate AF order persists up to 
larger doping concentrations in the electron-doped compounds as compared to the hole-doped ones 1 . Despite the fact 
that our calculations are done in the weak coupling limit, we believe that the main result still holds in the strong 
coupling case as well. In particular, the appearance of the Fermi surface pockets upon doping introduces a negative 
correction z to the spin stiffness. Its exact magnitude can be different in the strong coupling case; however, the 
asymmetry between the energy dispersions near the bottom of the conduction near (±71", 0) [(0, ±71")] points of the BZ 
and the top of the valence band near (±7r/2, ±7r/2) points of the BZ will be similar. Therefore, independent of the 
details we expect that the trend found in the weak-coupling, i. e. the relative stability of the commensurate AF order 
in the electron-doped cuprates as compared to the hole-doped ones, will remain. 

A natural question arises: when does the commensurate SDW order becomes unstable for the electron-doped case? 
Upon increase of the electron doping the pockets at (7r, 0) and (0, 7r) grow in size. However, the Pauli susceptibility 
in two dimensions is doping independent and consequently the doping does not immediately influence the commen- 
surability of the SDW order. At the same time, it is natural that for increasing x the SDW gap will decrease. This 
will reduce the gap between the a and the j3 band in principle, at certain doping level one would expect that both 
hole and electron pockets may appear at the Fermi level, as shown in FigJTJc). Although this situation seems to 
be supported by ARPES experiments in the electron-doped cuprates, the commensurate SDW order will become 
unstable in this case. For the same reason as in the hole-doped case, the presence of the hole pockets at (7r/2,7r/2) 
immediately introduces critical scattering for the commensurate spin waves. An example of this behavior is shown 
for the spin excitations in Figj3jc),(f) where the commensurate SDW order produces two types of the pockets, see 
FigjTJc). One immediately finds that spin stiffness is again negative and there are incommensurate peaks forming 
around Q. Therefore, the system will either develop incommensurate spiral SDW order or it will readjust the chemical 
potential such that the hole pockets will disappear with a corresponding change of the electron pocket's size. There 
is, however, another interesting possibility which we would like to discuss in the next section. 

II. SPIN EXCITATIONS IN THE COEXISTENCE AF+DSC PHASE 

As shown above, the presence of electron pockets at (±7r, 0) [(0, ±7r)] points of the BZ do not make the commen- 
surate AF order at (7r, 7r) unstable. This in principle opens an interesting possibility to have a nodeless d x i - y i -wave 
superconducting gap in the electron-doped cuprates which would coexist with AF order. Its features have been 
discussed previously in a number of mean-field studies^nEIl. Nevertheless, our results indicate that in many cases 
commensurate AF order becomes unstable with doping due to a negative spin stiffness introduced by the dopants 
carriers. Therefore, to see whether commensurate AF order can coexist with dSC requires a simultaneous analysis of 
the spin excitations. In addition, the study of the spin excitations itself in the unconventional superconducting state is 
interesting on its own due to formation of the spin resonance mode at the wave vector Q= (7r,7r) and Q = £l res . Such 
a feedback effect is quite interesting because its presence allows to confirm the phase structure of the superconducting 
order parameter even in the situation when thermodynamics reveals an isotropic s-wave behavior as in some of the 
iron-based superconductors^. In the coexistence phase of AF-fdSC, the situation is trickier due to the presence of the 
Goldstone mode at Q and the anisotropy of the spin fluctuations in the transverse and the longitudinal channel. 

To study the spin dynamics in the coexistence phase we follow ReffTSland apply the SDW Bogoliubov transformation 
to the pairing interaction [third term in Eq.dTl], and subsequently perform a mean- field decoupling in the particle- 
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particle channel, keeping anomalous expectation values of the form (a^ 



"*~ k |) and (/?kt^-k^)> anc ^ their complex 
conjugates. As in Ref. [18|we further neglect terms of the form (c^ 4,) due to the FS mismatch between the a 
and j3 bands. As we show below these terms can be safely ignored once the Umklapp Cooper-pairing terms are absent 
in the coexistence phase. The resulting mean-field Hamiltonian can be diagonalized by two independent Bogolyubov 

a, j3) as the energy dispersion of the two bands. The SC gaps 
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A£'^ are determined self-consist ently from two coupled gap equations, derived previousl}^ 
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where L k>p = (Vk-pf^ - V k _ p+Q F k ^) , M k>p = (v k _ p ^ - V k _ p+Q ^) with \ k /' k . 
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For the ^2 -wave pairing 



potential Vk-p = ^f^k^p case where 0k = cos/c x — cosk yj the SDW coherence factors can be factorized^ and one 
finds 
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Note, Dl is unity if \E^\ < fmj D , 
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bye frequency. The SC OP in the a and /3-bands are in phase without additional line nodes, in contrast to the s-wave 
cag ^i8H20] Nevertheless, the coexistence of commensurate AF and dSC phases can generate a higher harmonic, Ai, in 
the d x 2 _ y 2 -wave gap which has an opposite sign on the electron and hole FS pockets. This harmonic is proportional 
to the magnitude of the SDW and dSC gaps and arises due to Umklapp Cooper-pairing in terms of original fermions, 
i.e. averages of the type (ck,tC-k-Q,^)- These averages appear naturally in the coexistence phase as the wavevector 
Q becomes the new reciprocal wave vector of the lattice in the SDW state. At the same time, due to additional 
breaking of the spin rotational symmetry associated with the SDW transition the Umklapp Cooper-pairing terms 
formally belong now to the spin -triplet component of the Cooper-pair wave function with m z = as was discussed 
previously by several authors^ES. This indicates that the appearance of the Ai would be then associated with an 
additional phase transition in the dSC+SDW coexistence phase with a further change of the underlying symmetry of 
the mean-field Hamiltonian, i.e. from SU(2) x U(l) to 50(5) symmetry. 

For self-consistency the SC gap equations have to be complemented by the new equation which determines the 
SDW gap value in the dSC+SDW state 
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equations are complemented by the calculations of the total energy 
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to guarantee that their solution refers to its minimum. 

We further mention that a comparison of our eigenenergies Qj, with those found in Ref l22l bv the diagonalization of 
the 4x4 matrix allows us to make a connection between our approach of taking the sequential transformations and 
the exact diagonalization of the mean-field matrix. In particular, it is straightforward to see that the eigenenergies 
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FIG. 4: color online) Calculated mean-field phase diagram of coexisting <i-wave superconductivity and commensurate AF order 
in the electron-doped cuprates for the case Ai = 0. The solid and dashed blue curve refer to the Lifshitz transition indicating 
the change of the Fermi surface topology from two type of pockets to the one type of pockets in the pure AF and AF+dSC 
state, respectively. The open circle for x = refers to the AF semimetal with two equal sizes electron and hole pockets (high 
T) to AF insulator (low T) transition. 



are the same for Ai = 0. Thus, the use of sequential unitary transformations, performed by us, is exact as long as the 
triplet component of the Cooper-pair wave function is absent. For Ai ^ the eigenenergies agree only up to terms of 
the order ~ O(Ai) but start to differ for the higher order terms. To recover the same energy spectrum as in Refl22l 
for the coexistence phase one needs to take into account the Cooper-pair terms (a^ ^/?L k i ) in the AF state. They can 
be particularly important in the case when the magnetic gap becomes small when the interband Cooper-pairing may 
be non-negligible. Furthermore, the particular form of the pairing interaction in the momentum space can further 
modify the structure of the superconducting gap equations in the coexistence state. 

The expression for the transverse spin susceptibility in the coexistent state within a RPA remains structurally 
the same as in the pure AF state, see Eq.(3), except that the bare susceptibilities contain Cooper-pair creation and 
annihilation processes in addition to the quasiparticle scattering within one band and between the bands. The full 
expression is given in the Appendix. The situation is more complicated for the zz component as there is a mixing 
of the charge and longitudinal spin susceptibilities at Q for a finite dopingpl Nevertheless, the pole in total RPA 
susceptibility is still determined by 1 — Xo*(q, q, which allows us to use Eq.(7) also in the coexistence phase. 

In Fig(4] we present the results for the mean-field phase diagram of the electron-doped cuprates for coexisting 
commensurate AF and d-wave superconducting order. Although this type of phase diagram was already obtained in 
the literature, there are several features which are important to mention. Observe, for example that within a pure 
SDW phase at finite doping there is a Lifshitz transition (blue dashed curve) separating phases with a different FS 
topology with either one or two types of FS pockets. At higher temperatures, both electron and hole type of pockets 
are present at the FS, while below the Lifshitz transition only the electron pockets are present. Another interesting 
feature concerns the character of the phase transition into the coexisting AF+dSC state at low temperatures. In 
particular, analyzing the free energy we find that the transition from AF to AF+dSC state is of first order as a 
function of doping, while it becomes second order as a function of temperature. Furthermore, we notice that our total 
energy analysis shows that the stationary solutions for coexistence of SDW and d— wave superconductivity in the 
electron-doped cuprates have always slightly lower free energy in the case Ai = 0, i.e. when the triplet component 
of the Cooper-pairing is absent. In other words in our approach the symmetry of the problem remains SU(2) x U(l) 
in the coexistence regime. We believe that it is connected to the fact that we ignored the contribution from the 
interband Cooper-pair averages (^w/^Ik^) in the coexistence phase. Although small they could change the balance 
of the free energy towards the coexistence state with finite 'triplet' component of the Cooper-pairing. Furthermore, 
a modification of the momentum dependence of the Cooper-pairing interaction may also change the balance of the 
mean-field states. This, however, would require a separate analysis which is beyond the scope of the present paper. 

Next we compute the transverse spin response in the coexistence region of the phase diagram for three characteristic 
electron doping levels, n = 1.06, n = 1.09, and n = 1.12 (we fixed T= 0.02t). First it is easy to check analytically 
that as soon as A£ = A^ which is fulfilled for Ai = the condition for the pole in transverse susceptibility, 
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FIG. 5: color online) Calculated transverse, Imx^p A (q, q, ^) spin excitation spectra Q vs. q in units of tv/cl for three different 
electron dopings, n — 1.06, n — 1.09, and n — 1.12 (from upper to lower panel) for the coexisting SDW+dSC state and Ai = 
(right panel). For comparison the left panel shows the results for the pure SDW state. The blue lines denote 1 = £/RexJ _ (q, Q) 
condition. The intensity is shown on the log scale. The following parameters are used in the units of t for (a) \i — —0.5362, 
W = 0.5617, for (b) fi = -0.4573, W = 0.4617, for (c) fi = -0.3694, W = 0.3456, for (d) fi = -0.5349, W = 0.5530, 
A = 0.0727, for (e) /x = -0.4509, W = 0.4555, A = 0.0705, and for (f) fi = -0.3624, W = 0.3458, A = 0.0618. 



1 = L 7 "ReXn"~(q = Q, ^ = 0), again coincides with the mean- field equation determining the SDW gap in the coexistence 
phase, Eq.(ll). Therefore, similar to the pure SDW state the transverse component of the spin susceptibility shows 
gapless Goldstone mode also in the coexistence region. This is because the pure spin singlet d— wave superconductivity 
alone does not break the spin symmetry of the Hamiltonian. Here, our results do agree with those found very recently 
in a different formalism Ref,24, At the same time, in Refl25lthe spectrum of the transverse spin excitations computed 
numerically in the coexistence phase for hole doping seems to be gapped at the antiferromagnetic momentum, Q, 
which would contradict our analytical and numerical results. We speculate that this discrepancy arises from the fact 
that at the hole doping studied in Refl25lthe commensurate AF order is not a stable solution. In particular, although 
the mean-field solution for the commensurate AF order exists, the resulting spin waves have negative spin stiffness™ 
(see also our Fig. 3). As a result the true magnetic order in this case is incommensurate. 

Although the spin waves are gapless in the coexistence phase, their velocity is affected by the presence of d— wave 
superconductivity. In particular, in the FigJH] (d)-(f) we show the behavior of the transverse, Imx^p^(q, q, ft) spin 
excitations for three different electron dopings, n = 1.06, n = 1.09, and n = 1.12 (from upper to lower panel) in 
the coexistence state. Observe that in all three cases the spin excitations are gapless but the spin wave velocity is 
influenced by the superconducting gap. For small energies this can be seen analytically by expanding the denominator 
of Eqs. (4) around q,^ = up to quadratic order. This procedure leads to the spin wave velocity, c of the form: 



yt 2 (l/U-W 2 z) 
W 2 x 2 + (v)(l/U - W 2 z) 



(13) 



The coefficients for x : y, z contain in the lowest order terms proportional to Aq. Their exact form is given in 
the Appendix. The influence of the superconducting gap on the spin wave spectra is also seen by comparison of 
the transverse spin excitations spectra of the coexistence SDW+dSc state with that of the pure SDW state which 
is shown in the left panel, Figj5|a)-(c). The effect is particularly strong around 2Ao where the spin waves in the 
coexistence region exhibit a kink structure due to the interaction with the p — h continuum of the intraband a — a 
and j3 — j3 transitions. They are enhanced due to d— wave symmetry of the superconducting gap, i.e. due to the fact 
that one finds A^ = — A^ +q . for incommensurate momenta q^ > (0.8, 0.8)7r. This then leads to an enhancement of 
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FIG. 6: color online) Calculated longitudinal, Imx|fp A (q, q, Q) spin excitation spectra Q vs. q in units of ir/a for three different 
electron dopings, n — 1.06, n — 1.09, and n — 1.12 (from upper to lower panel) for the coexisting SDW+dSC state and Ai = 
(right panel). For comparison the left panel shows the behavior of the longitudinal susceptibility for the pure SDW state. The 
intensity is shown on the absolute scale. 



the intraband particle-hole continuum of both bands for Q w 2Ao- As the electron band around (±7r,0) and (0, ±7r) 
points of the BZ always crosses the Fermi level in the SDW state the enhancement of the particle-hole continuum of 
this band around 2Ao is responsible for the kink structure seen in the spin waves. In other words, the damping effects 
of the particle- hole continuum on the spin waves are present in both pure metallic SDW and coexisting SDW+dSC 
states. However, in the coexistence region there is also an effect of the strong renormalization of the spin wave due to 
the 2Ao structure of the particle-hole continuum of the intraband susceptibility, which then yields the renormalization 
of the spin wave velocity around 2Aq. Another interesting feature is that d— wave superconductivity stabilizes the 
commensurate AF state by partial gapping the particle-hole continuum in the coexistence state. Observe, for example 
that the spin waves computed for n = 1.12 in the pure SDW state show a tendency towards incommensurability, 
while in the coexistence state the spin excitations are still commensurate. 

In Figj6|d)-(f) we also show the behavior of the longitudinal response in the coexistence region. These are again 
compared to the pure SDW state Figj6^a)-(c). It is interesting that in this case the same structure of the particle- hole 
continuum at 2A for the intraband transitions which generates the strong damping of the spin waves around 2A in 
X +_ introduces the formation of the spin resonance at energies below 2Ao in \ zz - I n particular, the resonance occurs 
because of the discontinuous jump of the imaginary part of the intraband susceptibility at 2Aq. Correspondingly, 
there is a log singularity of the real part and below it Re xo(Q^) uo 2 for the intraband susceptibility. As a 
result, the total Re xo(Q?^) is enhanced at 2A and the resonance condition 1/U = Rexo(Q,^ = ^res) is fulfilled 
for Qres < 2Ao. This is in contrast to the pure SDW state where the main contribution at Q comes from the 
interband excitations which are gapped by the SDW gap magnitudes. At the same time, the intraband quasiparticle 
scattering in the SDW state is suppressed for the wavevector Q. In the coexistence state, in addition to the intraband 
quasiparticle scattering, there are also intraband Cooper-pairing annihilation processes, which start to contribute to 
the spin susceptibility at 2 A. These processes are enhanced by the d— wave symmetry of the gap and are responsible 
for the formation of the spin resonance. This resonance, however, now occurs in a very narrow q region, which 
is bounded by the dispersion of the particle-hole continuum. In particular, the spin resonance in the longitudinal 
component of the spin susceptibility is bounded from above at the wavevector Q by the onset of the particle-hole 
continuum associated with breaking of the Cooper-pairs at approximately 2Aq and from the left and right at finite £q 
by the onset of the continuum associated with intraband quasiparticle scattering which is gapless due to the presence 
of nodes on the a Fermi surface in the d— wave gap. As a result the spin resonance shows a characteristic upward 
dispersion and then becomes overdamped. This is quite different from the typical downward dispersion of the spin 
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resonance associated with the d-wave symmetry of the superconducting state without coexisting long range AF order. 

Finally we also note that our expression for the spin susceptibility computed for the metastable state for Ai ^ 
consistently signals its instability, that is we find 1 - UReXo~ (q = Q, ft = 0) < in this case. This indicates that the 
calculations of the spin response in presence of the triplet Cooper-pairing requires one to go beyond the sequential 
transformations, which we leave for the future studies. 

III. CONCLUSION 

In this work, we studied the behavior of the transverse and longitudinal spin excitation spectrum in a one-band model 
within a pure spin density wave (SDW) state and in the coexistence region of SDW and d-wave superconductivity. 
In particular, we analyze the evolution of the transverse and longitudinal spin excitations in doped AF metals. We 
have shown that for the sign of t' appropriate for the cuprates, commensurate spin excitations are stable within 
weak-coupling only on the electron-doped side of the phase diagram. Otherwise, the results we obtained in the SDW 
phase have been obtained in various forms by previous authors. However, in addition to the spin wave spectrum 
in this phase, we have focused on the evolution of particle hole excitations. These may play an important role in 
any instability of such an AF phase to superconductivity, since the pairing vertex due to exchange of the magnons 
themselves is suppressed due to the Adler principle. Calculations along these lines are in progress. 

Furthermore, motivated by the observed coexistence of the AF and ^2 -wave superconductivity in the electron- 
doped cuprates we computed the spin excitations in the coexistence region. We find that the Goldstone mode in the 
transverse channel remains gapless in the coexistence regime and that the excitations in the transverse channel are 
dominated by the renormalized spectrum of the spin waves. This happens, however, only if the Umklapp Cooper- 
pairing term, (ckC_k-Q), which actually belongs to the spin-triplet component of the Cooper- wave function is absent. 
We show also that within our approach this triplet component never appears as the free energy is always higher for the 
solution with non-zero Umklapp terms. This is a subtle point which requires a careful, fully self-consistent treatment 
of competing order parameters in the problem. At the same time, we find that the excitations in the longitudinal 
channel include a resonance mode at the commensurate momentum close to (7r,7r). The simultaneous coexistence of 
the longitudinal resonance and the transverse spin waves opens up an interesting possibility to use inelastic neutron 
scattering (INS) to identify the microscopic coexistence of the superconductivity and antiferromagnetism. We also 
stress the importance of well-defined particle-hole branches of the spin excitation spectrum which should also be 
observable in INS and should help to confirm our general picture. 

IV. ACKNOWLEDGEMENTS 

We acknowledge helpful discussions with G. Blumberg, A.V. Chubukov, R. Fernandes, J. -P. Ismer, R. Moessner, 
D.K. Morr, and J. Schmalian. W. Rowe and P.J.H. were supported by NSF-DMR- 1005625. W. Rowe is thankful for 
hospitality of the Ruhr-University Bochum where the final stage of this work was done. IE acknowledges financial 
support of the SFB Transregio 12, Merkur Foundation, and German Academic Exchange Service (DA AD PPP USA 
No. 50750339) and is thankful to the University of Florida in Gainesville where the project was initiated. JK 
acknowledges support from the Studienstiftung des deutschen Volkes and the IMPRS Dynamical Processes in Atoms, 
Molecules and Solids. 



* Electronic address: wwang@uil. edu| 

* Electronic address: ieremin@tp 3.rub.de| 

1 N. P. Armitage, P. Fournier, and R. L. Greene, Rev. Mod. Phys. 82, 2421 (2010). 

2 M. Inui, S. Doniach, P.J. Hirschfeld and A.E. Ruckenstein, Phys. Rev. B 37, 2320 (1988). 

3 B.I. Shraiman and E. D. Siggia, Phys. Rev. B 46, 8305 (1992). 

4 V. I. Belinicher, A. L. Chernyshev, and V. A. Shubin, Phys. Rev. B 56, 3381 (1997). 

5 A. Liischer, A.I. Milstein, and O.P Sushkov, Phys. Rev. B 75, 235120 (2007). 

6 P. J. Hirschfeld, M.M. Korshunov, I.I. Mazin, Rep. Prog. Phys. 74, 124508 (2011). 

7 C. Weber, K. Haule, and G. Kotliar, Nature Phys. 6, 574 (2010). 

8 J. R. Schrieffer, X. G. Wen, and S. C. Zhang, Phys. Rev. B 39, 11663 (1989). 

9 R. Coldea, S.M. Hayden, G. Aeppli, T.G. Perring, CD. Frost, T.E. Mason, S.-W. Cheong, and Z. Fisk, Phys. Rev. Lett. 
86, 5377 (2001). 

N.M.R. Peres and M.A.N. Araujo, Phys. Rev. B 65, 132404 (2002). 



12 



11 J. Knolle, I. Eremin, A. V. Chubukov, and R. Moessner, Phys. Rev. B 81, 140506(R) (2010). 

12 A.V. Chubukov, and D.M. Frenkel, Phys. Rev. B 46, 11884 (1992). 

13 Note that for t' ^ there is a solution for SDW gap at finite temperature such that it yields a compensated SDW metal 
with equal volume Fermi surface pockets of the electron and hole character. 

14 A.V. Chubukov, and K.A. Musaelian, Phys. Rev. B 51, 12605 (1995). 

15 A. Singh, Z. Tesanovic, and J.H. Kim, Phys. Rev. B 44, 7757 (1991). 

16 T. Das, R.S. Markiewicz, and A. Bansil, Phys. Rev. B 74, 020506(R) (2006). 

17 X.-Z. Yan, Q. Yuan, and C.S. Ting, Phys. Rev. B 74, 214521 (2006). 

18 J.-P. Ismer, I. Eremin, E. Rossi, D.K. Morr, and G. Blumberg, Phys. Rev. Lett. 105, 037003 (2010). 

19 L. L. Daemen and A. W. Overhauser Phys. Rev. B 39, 6431 (1989). 

20 L. N. Bulaevskii, A. I. Rusinov, M. L. Kulic, J. Low Temp. Phys. 39, 255(1980); M. L. Kulic, A. I. Lichtenstein, E. 
Goreatchkovski, and M. Mehring, Physica C 244, 185 (1995). 

21 G. C. Psaltakis and E. W. Fenton, J. Phys. C 16, 3913 (1983). 

22 B. Kyung, Phys. Rev. B 62, 9083 (2000). 

23 A. Aperis, G. Varelogiannis P.B. Littlewood B.D. Simons J. Phys.: Condens. Matter 20, 434235 (2008). 

24 H.-J. Lee, and T. Takimoto, rarXiv:1204.5032| (unpublished). 

25 T. Das, R. S. Markiewicz, A. Bansil, A. V. Balatsky, Phys. Rev. B 85, 224535 (2012). 



Appendix A: Spin susceptibility in the coexistent AF+dSC state 



Here, we present the expression for the spin susceptibility in the coexistent AF+dSC region. In particular the transverse 
component has the form 
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Here, + or — refers to the 7 = a and 7 = /3, respectively. Correspondingly, the longitudinal component has the form 
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To evaluate the spin wave velocity, we expand the denominator of Eqs. (4) around q,Q = up to quadratic order. This 
procedure leads to the spin wave velocity, c of the form 
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and the other one from the interband contributions 
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